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ABSTRACT 

Aims. In most microquasars, low-frequency quasi-periodic oscillations (LFQPO) have been classified into three types 
(A, B and C depending on the peak distribution in the PDS and the shape of the noise) but no explanation has 
been proposed yet. The accretion-ejection instability (AEI) was presented in 1999 as a possible explanation for the fast 
varying LFQPO that occur most often. Here we look at a possible generalization to explain the characteristics of the 
other two LFQPO types. 

Methods. It was recently shown that when the disk approaches its last stable orbit, the AEI is markedly affected by 
relativistic effects. We focus on the characteristics of the LFQPO that would result from the relativistic AEI and 
compare them with the different LFQPO types. 

Results. The effects of relativity on the AEI seem to be able to explain most of the characteristics of the three types of 
LFQPO within one formalism. 

Key words. X-rays: binaries, stars: individual (GRS 1915+105, XTE J1550-564, XTE J1859+226, GRO J1655-40, 
H17143-322), accretion disks 



1. Introduction 

Power density spectra (PDS) of black hole binaries show 
a high level of variability at all frequencies. These PDS 
are usually fitted with the sum of several Lorentzians. 
Depending on the spectral state of the black hole, thin fea- 
tures with a high power are also often detected. By a widely 
accepted convention these thin features are referred to as 
QPOs if the value of the coherence (Q — t'ccntroid/ FWHM) 
is higher than 2. Here we are interested in the LFQPO s 
with frequencies typically in the range 0.1 — 20Hz. These 
LFQPOs were originally considered to be one single phe- 
nomenon with varying frequencies. But during the 1998 
outburst of XTE J1550-564 , the LFQPOs displayed highly 
varying properties which led to classifying th em into three 
types labeled A, B and C (see for example Wijnands et 



Table 1. 



la me ±. Summary of the prop erties of the LFQPO types 
from |Remillard et al. (120021 ) and |Casella et q/.| ( |2Tj04| ), based 
on the microquasar XTE J1550-564. 



Properties 


type A 


type B 


type C 


frequency (Hz) 


~ 6 


~ 6 


0.1 - 20 


amplitude (%rms) 


3-4 


~ 4 


3-16 


Q v/FWHM 


~ 2 - 4 


rv, 4 


> 10 


phase lag (rad) 


-0.6 to -1.4 


to 0.4 


0.05 to -0.4 


subharmonic 




soft 


soft 


first harmonic 


soft 


soft 


hard 


coherence 


< 0.5 


~ 1 


~ 0.9 


HFQPO 


4/4 


6/9 


5/51 


noise 


weak red 


weak red 


strong flat top 



\aT. 1999 Remillard et al. 2002 1. Since then, these three 



types of LFQPOs have be en observed in seve ral sources, 
including GRS 1915+105 (|Soleri et al] [2007|, indicating 



that they are three genuinely different types of QPO, and 
that mechanisms common to all sources may give rise to 
them. 

Observational definitions of the three types of LFQPO 
are given m Tab le [Tj w hich is a summary of the results of |2002[ | Varniere et a/.[|2002[ ) 
Remillard et al.\ ( |2002[ ) and |Casella et at| ( |2004[ ) , based on 
tEe microquasars XTE J1550-564 and XTE J1859+226. 

The type C LFQPO is the one that occur most regularly. 
It is observed in a spectral state where the energy spectrum 
is close to a power-law with a photon index of about 1.5 to 



1.9 and with an exponential cutoff around lOOkeV, namely 
the low/hard state. Its main characteristic is a fast varying, 
highly coherent peak with a strong amplitude. It is observed 
simultaneously with a strong flat top noise in the PDS. It 
was proposed to be an expression of the accretion-ejection 
instability (AEI) i n conditions where it dominates the inner 
region of the disk (|Tagger fc Pellat 1999 Rodriguez et al. 



On the other hand, types B and A tend to appear in 
softer states such as the steep-power law state or soft in- 
termediate state, which shows blackbody and power-law 
components of similar amplitude. The main peak frequency 
of the PDS is not varying much compared to type C and 
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is generally observed in absence of a flat top noise in the 
PDS. Among their differences, the type B typically has 
a subharmonic as well as a first harmonic, while type 
A seems usually not accompanied by harmonics. Type B 
LFQPOs have a coherence similar to the type C, while type 
A are much less coherent. A recent study of the 'cathe- 
dral' LFQPO of XTE Jf859+226 (categorized as a type 
B) shows that the different peaks in the PDS might not 
be r elated to one another in just a simple harmonic man- 
ner ( Rodriguez fe Varniere l [2011 ). There is also a link 



between the type of LFQPO and the presence of a high- 
frequency QPO. In Table [1] we see that HFQPO are almost 
always detected when a type A or a type B LFQPO is also 
present, which is rarel y the case when type C LFQPOs are 
present. Furthermore, jRemillard et al. (2002) showed that 
when both LFQPO and HFQPO are present, their Q- values 
appear anti-correlated, which could hint at a competition 
mechanism between them. 

Up to now, theoretical models of LFQPO have focused 
on type C LFQPO and tende d to tie its frequency to a m ag- 



netoacoustical frequency (e.g. Titarchuk & Fiorito [2004 ) or 
to the Kcplcrian frequency at some radius in the disk (e.g. 



Tagger & Pellat 1999). The apparent different behavior of 
the three types of LFQPOs reposes the question whether 
they are coming from different mechanisms or if it is the 
same mechanism, only expressing differences occurring in 
the system. The first step in answering this question is to 
test if one mechanism can exhibit the different behaviors 
observed. Here we focus on the properties of the AEI and 
test if they could explain the peak distribution in the three 
types of LFQPOs in a single framework when relativistic ef- 
fects are taken into account. To do this, we choose to focus 
on an interpretation for the peaks and not for the contin- 
uum of the PDS, such as the band-limited noise (BLN) 
which sometimes occurs, or its correlations found with the 
main peak's frequency. We consider this correlation to be 
an expression of their origin and not an integral part of the 
QPO mechanism. 

In section [2] we briefly review the properties of the AEI 
that made us consider it as a good candidate for the most 
common LFQPO, the type C. We then investigate how rel- 
ativistic effects can modify the properties of the AEI when 
the inner edge of the disk approaches the last stable or- 
bit. The following section presents the LFQPO type that 
would result from the relativistic AEI (R-AEI). In the last 
section we turn to real data and see if the different flavors 
of the AEI can explain the observed characteristics of the 
three types of LFQPOs, first XTE J1550-564, in which the 
different types of LFQPO have been most frequently ob- 
served. Subsequently, we consider several other objects to 
gain an exhaustive view of the different types. 



where Q and k are the rotation and epicyclic frequencies 
(in a Keplerian disk Q = k), £ is the surface density and 
B is the equilibrium magnetic field. This criterion is ful- 
filled in disks with 'reasonable' density and magnetic field 
profiles. For the sake of completeness we present here the 
main prop erties of the AEI as they have been discussed 



elsewhere 
IVarniere et al 



Tag ger fc Pellat[ [ 1999 
2M3[ IVa rniere fc 



|Caunt fc Tagger |2001 
?agge r[ |2002| |Rodriguez 



et al. 2002 Varniere et all 120021 |Mikles et al. 2009) in 



association with the type C LFQPO - 

First of all, the AEI is an instability, therefore it can 
grow naturally to a high amplitude without a need for an 
external excitation. It is also able to account for the follow- 
ing observational characteristics: 

- The rotation frequency of the dominant m — 1 mode, 
i. e the one-armed spiral predicted by the AEI is a few tenths 
of the Keplerian frequency at the inner edge of the disk. 
This frequency is cons istent with the LFQPO frequency 



(Tagger & Pellat, 1999). 



Linear theory and nonlinear simulations 



Tagger 2001 ) show that the AEI forms a standing wave 



Caunt & 



pattern that can saturate at a finite amplitude. It can thus 
account for the persistence of a QPO whose quality factor 
should be limited only by the slow evolution of the factors 
(radial profiles of rotation, density, temperature, magnetic 
field) that fix its frequency and amplitude, or by nonlin- 
ear interaction between the modes when more than one is 
present. The result is therefore a pattern showing a fairly 
high coherence on time scales that are long compared to 
the frequency /orbital period. The AEI naturally explains 
the thin features (QPO) seen in the PDS, while broad band 
components (e.g. the band-limited noise) would be either 
QPOs broadened by other effects or totally different phe- 
nomena. 

- Bbased on variations in disk properties at the loca- 
tion of the spiral wave, the AEI partially reproduces the 
observed X-ray flux modulation ( Varniere et al.\ 2003). 

- The AEI can transfer to the corona a significant frac- 
tion of the energy and angular momentum that it extracts 
from the disk through Alfven waves, thus providing a sup- 
ply of energy that might feed the compact jet, often de- 
tected in the low-hard stat e, where the LFQPO is observed 
( |Varniere fc Tagger[ [2002| . 

- Taking into account the effects of General Relativity 
through the existence of a last stable orbit and orbital 
velocity profile, when the inner edge of the disk approaches 
its last stable orbit, we were able to explain the observed 
turnover in the correlation between the color radius 
(related to the inner disk radius, as d etermined by the 



spectral fits) and the LFQPO frequency dRodriguez et al. 
2002| |Varniere et al] [2002) |Mikles et al.\ |2rj0^| ~ 



2. The different flavors of the accretion-ejection 
instability 

2.1. The AEI as a model for the common varying LFQPO 



The AEI was first introduced by |Tagger fc Pellat | ( [1999 1 . 
It is a global instability occurring in disks threaded by a 
poloidal magnetic field on the order of equipartition with 
the gas pressure ('fully magnetized' disk), namely when the 
plasma /3 ~ 1 and it also requires 



d_ f K 2 S 
lh \2ilB 2 



>0, 



(1) 



2.2. The relativistic AEI 



In a recent paper, Varniere et al. (2011 hereafter paper I) 
showed that in a disk whose inner radius is close to the last 
stable orbit, the AEI could co-exist with another instabil- 
i ty, the Rossby wave insta bility (RWI), which was proposed 
( Tagger & Varniere 2006 ) as a possible explanation for the 
high-frequency QPO observed in microquasars. This would 
therefore create a state with both high-frequency and low- 
frequency QPOs, just as is observed for the type A and B 
LFQPOs in sources like XTE J1550-564. This motivated us 
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to investigat in more detail at the way the AEI is modified 
by relativistic effects, and how the observed flux modula- 
tions would be affected. 



of GRO J1655-40 during its 1996 outburst, while observa- 
tions of GRS 1915+105 regularly occupy both sides of the 
curve. 



2.2.1. Range of interest for the location of the inner disk edge 2.2.2. Numerical simulation results and mode spectrum 



For both instabilities (AEI and RWI) to coexist, we first 
need density and magnetic field profiles that meet the crite- 
ria Eq. 1. Additionally, the AEI requires the magnetic field 
and the gas pressure to be close to equipartition, namely 
[3 = 8irp/B 2 ~ 1, and the RWI requires the inner edge of 
the disk to be close to the last stable orbit, r in ~ Tlso- It 
was previously shown that the frequency of the AEI alone is 
modified by relati vistic effects when the disk approaches its 
last stable orbit (Varniere et aLl|2002|). This modification 



causes a turnover in the correlation between the frequency 
and the color radius, which is compatible with observations 



of GRO J1655-40 a nd GRS 1915+105 (see |Rodriguez et al. 
2002[|Mikles et a/.[|2009) . Fig. [I] shows the evolution of the 
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Fig. 1. Evolution of the frequency of the m — 1 mode of the 
AEI as a function of the position of the inner edge of the disk 
with respect to the last stable orbit. The RWI occurs when the 
inner edge of the disk is inside 1.3 tlso- See |Varniere et al.\ 
(2002) for more details. 

frequency of the m = 1 mode of the AEI as a function of 
the position of the inner edge of the disk with respect to the 
last stable orbit. We see that the range of ratio r in /rLSO m 
which both instabilities can coexist is the same in which the 
AEI would experience relativistic effects. In the right part 
of the curve in Fig. [I] the disk is mostly Keplerian, whereas 
in the left part relativistic effects near the inner disk edge 
modify the frequency-color radius correlation. Here we will 
focus on the top and left of the curve where we expect the 
relativistic AEI to have most of its impact. We have shown 



(Rodriguez et al. 2002 Mikles et al.\ 2009) that the left 



of the curve can correspond, for instance, to observations 



In paper I we presented numerical simulations of a disk 
subject to both instabilities (the RWI and the AEI) in 
its inner region, manifesting themselves as two spiral 
structures rotating at distinct frequencies. Although their 
results may depend in an unknown manner on the limits 
of the numerical model and on our hypotheses, they 
show quite stable features that we can compare with the 
observations. One of the main differences with simulations 
of the AEI alone is that instead of a dominating m = 1 
mode (where m is the azimuthal wavenumber, i.e. the 
number of arms of the spiral), we tend to have higher-m 
modes dominating for both the RWI and the AEI. Indeed, 
in the presence of the m — 2 mode of the RWI at the 
inner edge of the disk we tend to have a dominating m = 2 
mode of the AEI. These modes also tend to have a messier 
Fourier representation than for the AEI alone (see for 
example |Caunt &: Tagger) [2001 , for an example of the AEI 
alone). 

Because we are interested in the distribution of peaks 
in the PDS (to compare with LFQPO observations), we 
focus here on the evolution and strength of the m = 1 
to 4 contributions in the simulation, regardless of whether 
they are distinct modes or a mode and its harmonics. This 
does not allow us to compute the coherence of the differ- 
ent peaks, bur the AEI predicts quite a high coherence. 
Up to now, when trying to explain the LFQPO, the focus 
was on the strongest mode, here we deliberately looked at 
the evolution of the contribution up to the mode m = 4. 
Note that these modes are linearly independent. Their fre- 
quencies are eigenfrequencies of the system and depend on 
global properties of the disk. Thus linear theory shows that 
their frequencies are close to, but not exactly in, harmonic 
relations. In a PDS they would therefore appear as close to 
a fundamental and several (sub) harmonics. We are partic- 
ularly interested in their relative strengths, which can be 
directly compared with observation. 

Fig. [2] is a schematic view of the strength of the first 
four modes in the different configurations we found in the 
simulation. From left to right we see the evolution from a 
dominating m = 1 mode to a strong couple of (m = 2, m = 
4) modes. 

In Fig. [2^, we have a 'standard' distribution of the mode 
strength: the m = 1 mode is the strongest, followed by a 
weaker m = 2 mode, and even weaker m = 3 and m = 4 
modes. This would correspond to the peak distribution of 
type C QPO that can typically be seen in the low-hard 
state of GRS 1915+105 and also at the beginning of 
several outburst sources such as XTE J1550-564 and XTE 
J1859+226. 

Fig. [2)3 shows the case when the contribution of the 
m = 2 mode increases. We arrive at a balance between the 
strength of the m = 1 and the m — 2 modes; this would 
show a double peak in the PDS. In that case, the next 
mode in strength is the m = 4 mode and not the m — 3, 
as it was before. This corresponds to the specific c ase seen 
in XTE J1859+226 and dubbed cathedral QPO (iCasella 
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(a) 



(b) 



(c) 



(d) 



and ~ 9 Hz (see Mikles et al. 2009 for details), while 
for XTE J1550-564 the variation fortype B and A is ob - 
served to be between 5 Hz and 9 Hz (Remillard et al.\ 2002 1. 



Fig. 2. Shematic of the contribution behavior to the different 
mode of the AE1. (a) is a 'standard' peak distribution with 
higher modes of lower amplitude, whereas (b) is closer to the 
'cathedral' type. The distribution in (c) has been observed in 
type C and B LFQPOs, while (d) is more characteristic of type 
A LFQPOs. 



et al. 2004 Rodriguez & Varniere 2011) 



In Fig. [2}; the contribution of the m = 2 mode becomes 
dominant over the m = 1 mode. In a PDS this would 
appear as if the stronger peak harbored a weaker subhar- 
monic. The m — 3 mode is weak and of lower amplitude 
than the m — 4 mode but might be detected in bright 
sources. Depending on the PDS continuum and the shape 
of the noise, this would be labeled a type C (as in the later 
stage of outburst of XTE J1550-564 and XTE J1859+226.) 
or a type B once the source has reached a softer state. 

In Fig. [2ji this effect is even more accentuated because 
m = 4 becomes the second dominant mode after m = 2, 
while m = 1 is weak and m = 3 is barely detectable. 
This is similar to the type A QPO observed in outburst 
sources (see Fig. [4] for a fit of a type A LFQPO showing 
the two-peak structure). 



3. Characteristics of the R-AEI LFQPO and 
possible origin of the three types of LFQPOs 

From these simulation results we see that the X-ray signal 
from the disk can be expected to be modulated in distinct 
manners, well-suited to explain the observed differences 
between the three LFQPO types. 

- narrower frequency range : To show visible rela- 
tivistic effects on the AEI, the inner edge of the disk needs 
to lie close to the last stable orbit, which in turn limits 
the range of frequency as depicted in Fig. [T] Indeed, from 
linear analysis, it was shown that relativistic effects start 
to appe ar when the inner edg e of the disk is as close as 



3 r L so (Varniere et al. 2002[ ). When the inner disk edge 
is within 1.3 tlso the AEI exhibit stronger relativistic 
changes, as shown on the left of the Fig. [l] This limited 
radial range also limits the frequency range to a factor of 
a few, i.e. much less than in the nonrelativistic case where 
the inner edge of the disk can move much farther out, 
resulting in frequency variations by an order of magnitude. 
The narrower range of frequencies in this region is com- 
patible with the observed behavior of type A and type B 
QPOs, which we associate to the relativistic flavor of the 
AEI. Indeed, for GRS 1915+105 the observed frequency 
on the left of the curve in Fig. [I] varies between ~ 6 Hz 



- presence of a ' sub ' -harmonic : For the classical 
AEI, m = 1 is the dominant mode, which means that the 
lowest frequency also has the highest amplitude. The other 
observed contributions are much weaker, and given their 
weaker quality factor, will appear as a set of harmonics. 
For the R-AEI the dominant contribution is not m = 1 
but a higher m mode; for example the m = 2 mode of the 
AEI is the dominant contribution in presence of the m = 2 
mode of the RWI. In that case the weaker m = 1 mode 
could be interpreted as a subharmonic of the dominant 
772 = 2 mode. This would be a much more preferable 
interpretation because classically it is very difficult for a 
physical phenomenon to produce subharmonics. 

- appearance related with the occurrence of 
HFQP0 : Within the R-AEI framework we see that the 
LFQPO will be related to the HFQPO because both 
instabilities have similar conditions (see paper I for more 
details). 

Because of these characteristics the relativistic flavor 
of the AEI is a good candidate to explain the behavior 
of some type C as well as most of type B and A, while 
the nonrelativistic AEI can explain the regular type C (by 
regular we mean without a subharmonic) . There is no one- 
to-one relation between our model of mode evolution and 
the definition of the three types of LFQPO because we do 
not take into account the full PDS behavior, but only the 
LFQPO peak distribution. The phenomenological bound- 
ary between these three types takes into account other fac- 
tors such as the behavior of the continuum; these are not 
discussed here because they would require a full descrip- 
tion of the consequences of the instabilities on the disk, 
which is presently beyond our reach. At this stage, our in- 
terpretation thus only applies to and relies on the spectral 
distribution and evolution of the QPOs. 

type C LFQPO : Apparently, there are two cases for this 
type, as one can see in the top row of Fig. [3j The most 
common one is the evolving LFQPO without a sub- 
harmonic: this would be an expression of the AEI in 
its nonrelativistic form, i.e. the dominant contribution 
comes from the m = 1 mode as seen in Fig. [2]a,. In that 
state the inner radius of the disk is far from its last sta- 
ble orbit (larger than ~ 3 r^so)- But we also observe 
some type C LFQPOs, with subharmonics; we would 
explain them as an expression of the mildly relativistic- 
AEI with a dominant m — 2 mode. The main difference 
between this and type B resides in the underlying disk, 
namely a different spectral state. Therefore, the type C 
are represented as an evolution from case (a) to case (c) 
in Fig. § 

type B LFQPO : this would be an expression of a transition 
toward the relativistic AEI-dominated state with the 
m = 2 mode contribution dominating. Because the m — 
1 mode is present but not dominant, it would appear as 
a subharmonic of the dominant m = 2 mode. In that 
state the inner radius of the disk stays close to the last 
stable orbit (within 1.3 r^so)- In Fig. [2] type B would 
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0.1 1 10 0.1 1 10 

Frequency (Hz) Frequency (Hz) 




0.1 1 10 0.1 1 10 

Frequency (Hz) Frequency (Hz) 



Fig. 3. Top: two different configurations of the type C LFQPOs. Bottom: two different configuration of the type B LFQPOs 
depending on the relative strength of the two lowest peaks of the PDS taken during the 98 outburst of XTE Jf 550-564. 



be either case (b) or (c), as one can see from the PDS 
of XTE J1550-564 in the bottom row of Fig. [3j 
type A LFQPO : this would be an expression of the R-AEF 
dominated state with a strong contribution of both the 
m = 2 and m = 4 modes, while the modes m = 1 
and m = 3 are weaker. The case (d) of Fig. |] and the 
Fig. [4] are examples of this. Another important aspect 
of this case that we see in our simulations is that the 
modes are not as clearly defined as in the previous cases. 
This in turn would make the PDS appear 'messy' with 
a broader width for the detected QPO, as one can see 
from the Fig. [4] 



4. Comparison with observations 

We have shown that one instability was able to produce 
the main characteristics of the three types of LFQPO de- 
pending on the conditions in the accretion disk. We now 
investigate the observations in more detail, using several 
objects known to harbor different types of LFQPO to test 
our interpretation. 




Frequency (Hz) 



Fig. 4. Peak distribution of a type A LFQPOs taken during the 
f 998 outburst of XTE Jf550-564. 



4.1. The archetypical XTE J1550-564 

To compare the observational properties of the different 
types of LFQPOs with what is expected from the R-AE1, 
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we look first at XTE J1550-564, which exhibits all different 
types with regularity. In this case type A and B LFQPO 
have a small range of freque ncy (of about 5 Hz to 9Hz, 



see e.g. Remillard et al. 2002), which agrees with what we 



expect from the relativistic AEI 



Subharmonic or dominant "first harmonic" 

A striking feature of the LFQPOs during that outburst is 
the peak in the PDS at about half the frequency of the 
main peak. This feature was often noted as a 'subharmonic'. 
In the framework of the R-AEI presented here, we would 
explain this feature by the nondominant m — 1 mode. It 
is therefore important to determine if the dominant mode 
observed is possibly the m = 2 mode or if is indeed the 
fundamental. 





Fig. 6. Evolution of the rms amplitude of the two lowest peaks 
of the PDS. The circle represent the lowest peak at all times, 
while the square represents the peak detected at almost twice 
the lower peak frequency. The vertical line represents the first 
detection of an HFQPO in the PDS. Data taken from Remillard 
et fflZ.N2002l. 



Possible detection of the m = 3 mode 

To identify the first peak with the non-dominant m = 1 
mode, it would be interesting to look for a rather weak 
m = 3 m ode that the model predicts. Recently, |Rao et\ 



al. (2010 ) proposed that during this outburst up to four 



frequencies in the 1:2:3:4 ratio were seen, even though 
the third one is weak compared to the others. 

During several observations of this outburst, it was re- 
quired to add another peak to the fit at about three times 
the frequency of the first peak, as one can see in Fig. [7] and 
in the bottom right side of Fig. [3j Figure [7] is an observa- 



Fig. 5. Evolution of the frequencies of the two lowest peaks of 
the PDS o bserved during th e outb urst of XTE Jf 550-564. Data 
taken from |Remillard et al\ ( |2002[ >. 



To study the behavior of the two lowest peaks in the 
PDS, we took data from the 1998 outburst of XTE J 1550- 

((20021). Fig 



564 as they are published in Remillard et al 
[5] shows the evolution of the frequency of the two lower 
frequency peaks of the PDS as functions of time. The circles 
represent the lowest peak and the square the peak that is 
about twice (or more on a few occasions) the frequency of 
the lowest peak. The frequencies of the QPOs are, within 
the errors, compatible with a harmonic relation. 

In Fig [6] we see that both peaks are present throughout 
the outburst, even though their relative strengths change. 
This is coherent with an initially dominant contribution 
from the m = 1 mode, which is slowly replaced by the 
rn — 2 (and higher modes) as the disk approaches the last 
stable orbit and relativistic effects on the AEI becomes visi- 
ble. It is also interesting to note that around the time of the 
change of dominance for the different modes, the v — 183 
Hz HFQPO began to be detectable, reenforcing the idea 
that the disk is close to the last stable orbit and that rela- 



X 
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s 
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tivistic effects should play a role (Remillard et al. 2002). 



Frequency (Hz 



Fig. 7. Observation of XTE JI550-564 on MJD 51070.13 one 
can see frequencies near the 1:2:3 ratio. 

tion on MJD 5107.13, i.e. at the beginning of the outburst. 
At that date the strongest peak is also the one with the 
lowest frequency and we detect up to ?>v in the PDS. The 
fit requires seven Lorentzians (x 2 — 1.2 for 248 dof), three 
of which are narrow features and can truly be considered 
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as QPOs, and gives the following parameter for the peaks: 
Peak 1 = 3.185 ± 0.007 Hz, rms = 13.1±° A 3 %), 
Peak 2 (v = 6.32 ± 0.02 Hz, rms = 5.5±g;|%) and 

Peak 3 (u = 9.29±£;£g Hz, rms = 1.1 ± 0.3%). 

Thus for that observation we have significant peaks in the 

1:2:3 ratio in the PDS. 




Frequency (Hz) 

Fig. 8. Later during the outburst of XTE J1550-564 on MJD 
51099.6 we see 1:2:4 and a possible feature where the 3 would 
be. The arrow indicates the position of this potential thin 
feature which is not accounted for in our fit. 

On several other dates during that outburst a feature 
in the PDS around the m = 3 frequency seems possible, 
but is not required to obtain a good fit. Fig. [8] shows such 
an example; a thin feature, indicated by the arrow, may 
be present, although it is not statistically required by the 
fit. Should this feature be a genuine one, it would then be 
coherent with the presence of a weak m = 3 mode that is 
clearly detectable when the m = 1 mode dominates and 
becomes a weak feature as the m = 2 mode becomes dom- 
inant. 

Link with the mode of the HFQPO 

As was shown in paper I, the R-AEI happens simultane- 
ously with the the RWI, which was proposed to be at the 
origin of the HFQPO. This implies a strong link between 
the occurrence of type A and B LFQPO and the HFQPO. 
As we showed in Fig. [6j the dominant peak changed from 
the lowest frequency to the first harmonic on the same day 
that the HF QPO occured. It was also noted in Rcmillard 
et al. (2002) that the data suggest l that there is an anticor- 
relation between the amplitudes and coherence of HFQPOs 
and LFQPOs, respectively.' This point is in favor of a link, 
and possibly a competition, between the mechanism at the 
origin of both the HF and LFQPOs of type B and A such 
as the one presented here. 

When looking at the statistics of the occurrence of 
HFQPO with the different types of LFQPO, we see that 
the HFQPOs have frequencies near 184 Hz and 276 Hz, 
which appear while we have the type B (in six out of the 
nine observations of type B there was an HFQPO around 
180Hz) and type A (four out of the four observations of 



type A there was an HFQPO around 280Hz) of LFQPO 
respectively. 

Following our interpretation according to which the 
HFQPO is due to the RWI and the LFQPO is due to 
the AEI, we expect to sometimes see the to = 1 mode of 
the RWI, especially during the 'transition' from the mode 
rn = 1 to the mode m = 2 of the AEI. This can be expected 
since the strength of cither mode depends on many disk 
properties that affect their growth rates, so that the usual 
dominance of the m=2 cannot be an abo slute rule. During 
the 1998 outburst, [Remillard et al.\ j2002| detected at least 
one occurrence of an HFQPO at 92 Hz with a significance of 
3.2cr and an amplitude of 0.64±0.10%. A peak was detected 
at that frequency in several other observations but its Q 
value did not always qualify it as a QPO. The frequency of 
this peak is coherent with it being the 'to = 1' associated 
with the 2:3 ratio observed in that source. This could be 
the to = 1 mode of the RWI instability associated with an 
LFQPO where the relativistic effects are becoming signif- 
icant. Some more detailed study of this potential to = 1 
mode of the HFQPOs is required to properly compare it 
with the model. But this is still difficult at the moment 
because we have so few occurrences of this. 



4.2. XTE J1859+226 

The outbu r st of XTE J1859+226 in 1999 was studied by 
Cui et al.\ ( |2000 [), who found HFQPOs. They also noted 
a subharmonic 



to the LFQPO simultaneously with the 
HFQPO. Thi s hinted towar d diffe rent types of LFQPOs. 
Subsequently, Casella et al. (2004) studied the same out- 
burst, focusing on the timing. The detection of the three 
types of LFQPOs was confirmed and showed similar be- 
havior to the one observed in XTE J1550-564. Nevertheless, 
Casella et al. observed some interesting facts that were not 
seen in XTE J1550-564. 



The lowest peak stays the same 



Casella et al. ( 2004 1 noted that during this outburst a sub- 
harmonic was always detected except for the first observa- 
tions. Interestingly, the frequency of the only peak detected 
in the first observation is coherent with the frequency of the 
lowest peak (which they called 'subharmonic') of the sec- 
ond observations. This would suggest, in our model, that 
the m — 2 mode became dominant early in the outburst. 

More importantly, it was also noted that for the dura- 
tion of the outburst the lag related to the lowest peak is 
always negative, with a behavior independent of the type 
of LFQPO. This suggests a common origin for the lowest 
peak in the PDS, independently of the LFQPO type. In our 
model, it would correspond to the to — 1 mode of the AEI. 

A possible feature in 1:3 ratio with the lowest peak 



This outburst also exhibits case (b) of Fig. § that |Casella 
et al.\ ( |2004[ ) called 'cathedral' LFQPO (see Fig. [9) , where 
the two lowest peaks have similar amplitudes. A complete 
timing analysis of the two observations showing this pe - 
culiar type is presented in Rodriguez & Varniere (2011). 
Here we use Obs. 40124-01-24-00 to illustrate our model. 
The broadband 2.55-40 keV PDS was fitted with a sum of 
three broad and three narrow Lorentzians. A constant (in 
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standard Leahy normalized units) was added to account 
for the white noise. This is represented by the linear line 
in Fig. [9] where the plot is in power*frequency units. The 
three thin features have the following parameters: 
Peak f : F= 2.94 ± 0.02 Hz, Q=5.9 , A= 2.8 ± 0.1% 
Peak 2: F= 5.828 ± 0.025 Hz, Q= 7.3 ; A=4.7 ± 0.1% 
Peak 3: F= 11.2±g j Hz, Q= 9.5 , A=l.l ± 0.3 % . 
The three peaks therefore have frequencies that are corn- 
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Fig. 9. Twin peaks of the cathedral subclass of LFQPO type B 
in XTE J1859+226 with an harmonic around 11Hz and a broad 



feature around 9Hz (Rodriguez & Varniere 2011 1 



patible (at the ~ 2er level) with them being harmonics in 



the ratio 1:2:4 or 1/2:1:2 (Rodriguez k. Varniere 2011 1. 



Interestingly, there is an additional feature in the PDS at 
the expected frequency for an harmonic relation with peak 
1 in a ratio 1 : 3, namely F=9.0°q 5 . This feature is broad 
(Q=1.5) but close to the QPO limit and has an amplitude 
of A=3.2±0.2%. It is therefore possible that we are indeed 
observing the m = 3 mode along with the m = 1,2, 4. XTE 
J1859+226 is therefore the second object (after XTE J1550- 
564) to show a feature in which the weak m = 3 mode 
should be detected. Detecting this possible feature in an- 
other object strengthens its potential association with a 
weak r7i = 3 mode, therefore putting the lowest peak in the 
PDS as the 'fundamental' and not a subharmonic. 



Transition and competition between the types of LFQPOs 



During that same outburst, Casella et al. (2004) observed 
for the first time several transitions between the different 
types of LFQPO within the same observation. They no- 
ticed the absence of a direct transition from type A to type 
C LFQPO, which seems to emphasize the central role of 
type B LFQPO in the transitions. At the same time, they 
also noted the apparent 'balance' between the amplitude of 
the different peaks of the PDS. This is particularly strong 
regarding the compared strength of the 'subharmonic' and 
the "first harmonic". Indeed, at the strongest rms ampli- 
tude of the "subharmonic" the 'first harmonic' is barely 
detectable. The reverse is also true. This could be an ex- 
pression of some kind of competition between the two. 



This agrees well with our model, in which we associate 
each peak with a mode of the instability. This type of 
'balance' is similar to the evolution seen between (c) and 
(d) of Fig [2j Depending on the condition in the disk, 
the dominant mode will change and therefore create this 
transition. 



4.3. Other objects in which type B can be inferred 

Using observations from XTE J1550-564 and XTE 
J1859+226, we were able to find a good agreement between 
the characteristics of the R-AEI and the observations of 
type A and B LFQPO. We now checked if we can infer the 
presence of these types of LFQPO from the other charac- 
teristics shown here such as HFQPO or the LFQPO lying 
in the left part of Fig [TJ 

GRO J1655-40 

The LFQPO data of the 1996-97 outburst of GRO J1655-40 
were not classified into the three types as is done for more 
recent observations. Nevertheless, this object seems a good 
candidate to harbor type B and A LFQPO based on sev- 
eral facts. First of all, not only HFQPOs were discovered 
in the source, but the LFQPO and the spectral parameters 
have a peculiar correlation when compared with other mi- 



croquasars such as XTE J1550-564 (Sobczak et al. 2000). 
Indeed, the correlation between the LFQPO frequency and 
the color radius is inverted compared to other sources. This 
behavior was later explained by using the AEI as a model 
for the LFQPO and requiring that the inner radius of the 
disk was close to its last stable orbit (|Rodriguez et al. 2002 



Varniere et al.\ |2002 ) during most of the outburst. Indeed, 
GRO J1655-40 was the first source to be discovered on the 
left side, namely r int < l.Sr^goi of Fig. [T] This makes it, 
in our interpretation, a manifestation of the R-AEI in the 
disk. 

Looking more carefully at this outburst, we see that 
the LFQPO exhibits several of the characteristics of type 
B LFQPO. Indeed, it is shown to have an almost first- 
harmonic ('almost'because the two frequencies are not ex- 
actly in an harmonic relationship but very close to it) and 
sometimes that h armonic had a higher rms amplitude than 



the fundamental (ISobczak et all 2000). All of this make us 



consider GRO J 1655-40 as another candidate in our search 
for the different types of LFQPOs. Indeed, we looked at 
several dates showing those characteristics (HFQPOs, two 
peaks in close harmonic relationship) and confirm the pres- 
ence of non-C LFQPOs in that source. 

GRS 1915+105 

A high-frequency QPO was discovered in 
GRS 1915+1 05 during the 9 class of variab i lity by |Belloni| 



et al. (2006); more recently |Mikles et al (2009) showed 
that during the same class the color radius-frequency 
relation was different from the standard Keplerian one, 
and that the points were on the left side of the theoretical 
plot of Fig. [I] It would be interesting to study the LFQPO 
in more detail during the part of the observation where the 
LFQPO lies on the left of the correlation and see of what 
type they are. However, this is difficult because the source 
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is highly variable and spends little time on the left side of 
the theoretical plot of Fig. [T] 

The search was made in other classes of the source by 
Soleri et al. (2007), who found a transient type B LFQPO 
in the classes f3 and fi. It would be interesting to look at the 
9 class where the source spends more time in the top/left 
side of the Fig. [T] 

4.4. H 1743-322 and the search for HFQPO 

In the model presented here, the presence of type A 
and B LFQPOs happens in conjunction with HFQPOs. 
Observations tend to show a link. It therefore seems possi- 
ble to use the type of LFQPO as a tracer to narrow down 
the search for HFQPO in microquasars. 

This idea has already been used observationally. Indeed, 
in 2005 Homan & Belloni discovered a pair of HFQPOs in 
H1743-322 in a possible 3:2 ratio and mentioned the pos- 
sible presence of the '1' associated with the 3:2 HFQPOs. 
They also noted that, in presence of those HFQPOs, the 
LFQPOs where different from the 'standard' C type and, 
from their frequency, shape and harmonic content, found 
them to be closer to the A or B type as defined in |Wijnands| 
et al. ( 1999 ). They also reported in one of their observations 
(labeled obs 2. in their paper) what they interpreted as a 
transition from a type B to A LFQPO, which would link 
the two types of LFQPO more closely. 

Subsequently, |Remillard et al\ ( |2006[ ) tried to detect 
more occurrences of these HFQPOs. Following the link that 
seems to exist between HFQPO and the type of LFQPO, 
they decided to focus the search on observations that ex- 
hibit type B and A LFQPOs, namely LFQPO without 
band-limited noise. To test the link between the different 
types of LFQPO and the presence of HFQPO, they used 
the shape of the PDS to sort out the 130 observations of 
that source in four groups. The first group has no LFQPO, 
the second group exhibits 'standard' C type LFQPO with 
a band-limited-noise, and the two last groups had a non-C 
type LFQPO, the distinction between the two last groups 
being made depending on the flux level. They only found 
HFQPOs in these two latter states, confirming a link be- 
tween these non-type-C LFQPOs and HFQPOs. 



5. Conclusions 

We have presented the relativistic flavor of the AEI as a 
possible cause for some of the distinct characteristics for 
the type A and B LFQPOs in microquasars, with some of 
the semi-relativistic effects being visible already in some 
type Cs. The AEI evolves from its Keplerian to its rela- 
tivistic form when the fully magnetized disk extends close 
to its last stable orbit. The effect is strongest when the 
disk extends down to its last stable orbit. In these condi- 
tions two instabilities can occur, the AEI in its relativistic 
form, at the origin of the LFQPO, and the RWI, at the 
origin of the HFQPO. The relativistic flavor of the AEI is 
a good candidate to explain the behavior of type A and B 
of the LFQPO such as their relatively small amplitude in 
frequency, their (sub)harmonic content or the relationship 
with the HFQPO. 

It can also be seen as a cue in the search for the elu- 
sive HFQPO. Indeed, the behavior of the LFQPO, which is 
much easier to detect, can be taken as an indicator of the 



presence of HFQPO in the disk. In a first step, we were able 
to confirm that association by considering observations of 
GRO J1655-40 during its previous outburst beforereturning 
to other objects. 
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